Viscous corrections to the resistance of nano-junctions: 
a dispersion relation approach 
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It is well known that the viscosity of a homogeneous electron liquid diverges in the limits of zero 
frequency and zero temperature. A nanojunction breaks translational invariance and necessarily 
cuts off this divergence. However, the estimate of the ensuing viscosity is far from trivial. Here, 
we propose an approach based on a Kramers-Kronig dispersion relation, which connects the zero- 
frequency viscosity, 17(0) , to the high-frequency shear modulus, /j.^, of the electron liquid via 77(0) = 
with r the junction-specific momentum relaxation time. By making use of a simple formula 
derived from time-dependent current-density functional theory we then estimate the many-body 
contributions to the resistance for an integrable junction potential and find that these viscous effects 
may be much larger than previously suggested for junctions of low conductance. 
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I. INTRODUCTION 



Viscosity, namely the effect of momentum transfer be- 
tween adjacent layers of a liquid, is a fundamental con- 
cept in both classical and quantum physics [H-Q- In the 
case of the electron liquid, it was shown by Abrikosov and 
Khalatnikov (AK) more than fifty years ago that the vis- 
cosity of a homogeneous liquid diverges at zero frequency 
and zero temperature (in this precise order of limits) Q . 
The physical reason for this divergence is related to the 
fact that at zero temperature the quasi-particles close 
to the Fermi energy are very long-lived and can trans- 
port momentum to distances arbitrarily far away from 
the source of the perturbation. We thus expect that any 
potential that breaks translational invariance would cut 
off such divergence giving rise to a finite (albeit possibly 
large) d.c. viscosity at zero temperature. 

The problem has garnered renewed attention in the 
past few years in the context of transport in nanoscale 
systems [5{ , where it was shown explicitly that the Lan- 
dauer formula for the single-particle resistance of a nano- 
junction R s — h/(2e 2 J^j Tj) (the sum is over the eigen- 
channels of transmission Tj) @i f ans to include cer- 
tain many-body effects, which cannot be described as 
single-particle scattering from an effective potential. 0,0 
Within the framework of the time-dependent current 
density functional theory [j| Hoj . one clearly sees that 
such effects arise from the frequency dependence of the 
exchange-correlation (xc) field. [ll[ More precisely, one 
can split the full exchange-correlation potential into a 
static component, which controls the transmission proba- 
bility Tj of the Landauer formula, and a dynamical com- 
ponent, which corrects the Landauer formula In- 
deed, even though the control of nanoscale junctions at 
the atomic level is far from being ideal, it is now clear 
that the Landauer formula computed within the frame- 
work of ground-state density- functional theory (DFT), 
overestimates the measured conductance by at least an 
order of magnitude in the case of low-conductance struc- 



tures (e.g., molecular structures), while it provides rea- 
sonable agreement in the case of metallic quantum point 
contacts, which show high-conductance values |12j . Part 
of this discrepancy has been attributed to errors in de- 
termining the position of the energy levels of the system 
relative to the electro-chemical potential in the leads - 
errors which in turn are related to self-interaction cor- 
rections, discontinuities in the xc potential as a function 
of particle number, and so on. (l3j However, these correc- 
tions do not fully solve the discrepancy between theory 
and experiments: other many-body effects, in particular 
those related to the viscous nature of the electron liquid, 
may play an important role. 



II. DYNAMICAL CORRECTIONS TO THE 
RESISTANCE 

Time-dependent current density functional theory pro- 
vides some insight into the physical character of these 
many-body corrections. Indeed, it was shown that the 
dynamical corrections to the xc field give rise to a vis- 
cous force [9|, llOj , similar to the ordinary hydrodynamic 
viscous force, but controlled entirely by electron-electron 
interactions - the ordinary hydrodynamic viscosity does 
not explicitly depend on interactions, but relies implicitly 
on the presence of collisions capable of establishing a lo- 
cal thermodynamic equilibrium [1]. Thus, the resistance 
in excess of the Landauer formula could be interpreted 
as the effect of the extra "friction" arising from the xc 
viscosity. In Ref. Q, this additional many-body resis- 
tance - which was termed "dynamical" precisely because 
it vanishes in a strictly ground-state formulation of the 
theory [111 ] - was estimated (assuming no current den- 
sity variation in the junction) in terms of the exchange- 
correlation (xc) viscosity of the liquid @, Q 
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where x is the direction of current flow through the junc- 
tion, _L represents the transverse directions, S is the 
cross-section area of the nano-junction, n is the ground- 
state electron density calculated, e.g., from the self- 
consistent static density functional theory (DFT) with 
the xc functional treated within the local-density ap- 
proximation (LDA), d x and d± are derivatives along the 
current flow and perpendicular to it, respectively. Here, 
Vxc{0) > is the zero- frequency xc viscosity of the ho- 
mogeneous electron liquid at the average electron density 
of the junction, and e is the charge of the electron. Thus 
the dynamical correction increases the total resistance of 
nano-junctions to R to t = R s + Rdyn [HI- Notice that 
the non-uniformity of the electron density is essential to 
the effect: i?dyn vanishes, as it should, if n is spatially 
uniform. 

The precise definition of r) xc (0) in current density func- 
tional theory is 16] 

Vxc{0) = -n hm , (2.2) 

u— >0 UJ 

where f xc .T{^>) is the transverse component of the dy- 
namical exchange-correlation kernel of the homogeneous 
electron liquid at density n. The kernel itself is defined 
as the difference between the inverse current-current re- 
sponse functions of the interacting and non-interacting 
system at the same density. This quantity admits a per- 
turbative expansion in the interaction parameter r s of 
the electron liquid - the average distance between elec- 
trons expressed in units of the Bohr radius. Our initial 
estimates of i?dyn were based on an extrapolated high- 
density expansion of rj xc , which gave Rd yn of about 10% 
of R s for molecular junctions, but considerably smaller 
for metallic quantum point contacts @, Since then, 
the use of a more accurate expression for r\ xc as a function 
of r s , namely [HI ] 

hn 

60r s A/1 + 807 1 - 40r s 2/i + 62r s 1/6 

has been shown to produce considerably smaller correc- 
tions for the case of two infinite jellium electrodes sep- 
arated by a vacuum gap. [l7j Indeed, this r\ xc is of the 
order of 10~ 7 Joule-sec/m 3 for typical metallic densities 
such as r s — 3 for gold, thus naively suggesting that these 
dynamical corrections are small under all circumstances. 

Before jumping to conclusions, however, it must be 
noted that the expression (|2.3|) is not the appropriate 
zero- frequency viscosity, as required by Eq. (I2.2[) . but 
rather was calculated in Ref. [16j under the implicit as- 
sumption uj 1/r where r is the momentum relaxation 
time for a quasi-particle. This condition is easily sat- 
isfied in the homogeneous electron liquid in the limit of 
zero temperature T, since r (in the absence of impurities) 
tends to infinity as 1/T 2 . It is certainly not satisfied in a 
nano junction where the main limiting factor to the quasi- 
particle lifetime is the geometric size of the constriction 
itself. [!, [Tsl - [20l j Indeed, it is precisely the geometrical 



constriction experienced by the electron wave packets as 
they move into the nanojunction that introduces a short 
momentum relaxation time 0, [H| , which in turn cuts off 
the divergence of the viscosity of the uniform electron 
liquid. 

We are then led to consider the opposite and physi- 
cally correct limit of uj <C 1/r. But here, we run into 
the problem that the non-uniformity of the electron liq- 
uid must be fully taken into account. If, for instance, 
one calculates the d.c. viscosity of the homogeneous elec- 
tron liquid, using the techniques developed by AK Q one 
obtains (employing the Thomas- Fermi approximation for 
the screened electron-electron (e-e) interaction) 

hn ,1.813 x 10 3 , 2 r tt(1 + 2r s ) ?n „ s 

^ = (^ ( ^— Hariri} ,(2 - 4) 

where f s = ar s /n and a = (4/97T) 1 / 3 ~ 0.521. This 
expression has two major shortcomings. First and fore- 
most, it diverges at low temperature as 1/T 2 , which 
is the consequence of undisturbed momentum trans- 
port by the long-lived quasi-particles of the homogeneous 
electron liquid to distances arbitrarily far away from 
the source of perturbation. Even at room temperature 
Vak = 1.12 x 10 4 /T 2 Joulc-sec/m 3 (temperature T is 
in Kelvin) for gold (r s — 3), which, if used in Eq. (|2.1[) 
would produce unreasonably large values of i?dyn- Sec- 
ond, the AK viscosity is conceptually different from the 
xc viscosity: it arises from the full current-current re- 
sponse function, not just from the exchange-correlation 
kernel, and there is no simple way to separate the latter, 
since the AK calculation is non-perturbative with respect 
to the electron-electron coupling strength. Instead, in our 
case, such a separation is essential, since the larger (elas- 
tic) part of the resistance is already taken into account 
exactly by the Landauer formula, and only the xc ker- 
nel contributes to the dynamical (inelastic) correction. 
Thus, we need a way to "zero-in" on the xc viscosity (de- 
fined formally in terms of the xc kernel f xc in Eq. (|2.2[1 ) 
just as we did in Ref. [l6| - but now we must go to the 
opposite regime of uj <C 1/t. 

Equipped with this understanding, we first show that 
the viscosity of the electron liquid in the nanojunc- 
tion can be estimated by a dispersion relation approach, 
which combines information about the high-frequency 
elastic properties of the interacting electron liquid with 
the short momentum relaxation time induced by the con- 
striction. We then estimate the dynamical contribution 
to the resistance of nanostructures from Eq. (|2.1[) and 
show that this contribution is relatively small for nearly 
transparent junctions but becomes sizeable in the limit 
of zero transmission. 

III. DISPERSION RELATION APPROACH 

Let us then start by recalling that the viscosity tj(uj) 
and the shear modulus /i(w), regarded as functions of fre- 
quency, are respectively the imaginary and the real part 
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FIG. 1: The comparison between the viscosity calculated by 
Abrikosov and Khalatnikov, tjak, and the viscosity calculated 
using the dispersion relation (|3.2p for different values of r s . 
The connecting lines are a guide to the eye. 



of a visco-elastic modulus fi(uj) = fi(uj) — iu)i](uj). [Il.l3l.llll] 
The shear modulus vanishes at zero frequency, since we 
are in a liquid state, and tends to a finite value fj,oo a t 
infinite frequency, where by "infinite" we mean a fre- 
quency much larger than 1/r. On the other hand, the 
viscosity, rj(uj), varies from the desired value, 77(0), at 
zero frequency to zero at "infinite" frequency. Then, the 
Kramers-Kronig dispersion relation tells us that 



= jLt(O) = /loo / Tj{uS)<ko . 

I" Jo 



(3-1) 



We do not know the detailed frequency dependence of 
T}{uS), but we can assume that it roughly switches from 
the d.c. value 77(0) for oj < 1/r to approximately zero for 
uj > 1/t. This simple reasoning leads from Eq. (|3.ip to 
the relation d !(} 



?7(0) ~ fi a 



(3.2) 



The importance of this relation, in the present context, 
is that it connects the quantity of primary interest, the 
zero-frequency viscosity, to two quantities that can be 
rather easily estimated, namely, the high-frequency shear 
modulus - a positive definite quantity that can be ex- 
pressed in terms of the exact energy of the electron liq- 
uid - and the momentum relaxation time r, which, as 
we have discussed above, in a nano-junction is mostly 
controlled by elastic boundary scattering from the con- 
fining geometry. A crucial feature of this relation is 
that it allows a neat separation of 77(0) into two parts: 
a "single-particle" contribution, which arises from the 
non-interacting kinetic part of the shear modulus (see 
Eq. (|3.3p below), and an "exchange-correlation" contri- 
bution, which actually determines, via Eq. (|2.ip . the dy- 
namical correction to the resistivity. 

As a first test of the reasonableness of this approach 
and also to lend support to its quantitative accuracy, let 
us apply it to the calculation of 77(0) for the homogeneous 



electron liquid, with fi^ the full high-frequency shear 
modulus given, for example, in Eq. (103) of Ref. [2lj ]. 
This can be written as 



fJ-s 
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f^xc 
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(r s ) 



(3.3) 



where fi s = 2ne^/5 is the non-interacting shear modulus 
(arising from Pauli exclusion principle) with tp being 
the Fermi energy, and the remainder, /i X c,<x>, is expressed 
in terms of the exchange-correlation energy per particle, 
e xc , and its derivative eL c with respect to the electron 
liquid parameter r s . [22| In the homogeneous electron 
liquid the momentum relaxation time, r, is limited only 
by electron-electron interactions and therefore diverges 
with temperature T as 1/T 2 :0] 



1 
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&(r.) (3.4) 



where £k is the energy of the quasi-particle. For a typical 
density of r s = 3, £ 3 (r s = 3) ~ 0.5 and r ee ~ 5 x 10~ 7 /T 2 
s at the Fermi energy. For this set of parameters, from 
Eq. we then obtain rj(0) = ~ 1.6 x 10 4 /T 2 

Joule-s/m 3 which is in good agreement with tjak = 
1.12 x 10 4 /T 2 Joule-s/m 3 . The same level of agreement is 
found for a wide range of r s values of typical metallic sys- 
tems, as shown in Fig. Q] It is also worth stressing that 
the agreement we find with our dispersion relation ap- 
proach is extremely good for r s < 4 despite the assump- 
tions we have made in reaching Eq. (|3.2p . The reason why 
the agreement between the two viscosities decreases with 
increasing r s can be attributed to the fact that Eq. (13. 4p 
for the relaxation time is strictly valid for small values of 
r s . These results thus lend strong support to the present 
approach to compute viscosity in terms of the relaxation 
time t. 



IV. AN INTEGRABLE MODEL 



Having checked that the dispersion relation (|3 . 2[) gives 
both a qualitative and a quantitative account of the d.c. 
viscosity we can now proceed to the estimate of the lat- 
ter in the presence of a nanojunction. We consider the 
following confining potential (see Fig. [2] for its schematic) 
which - despite its simplicity - is a reasonable choice to 
mimic the scattering properties of nanostructures 



V(x,y,z) 



A + B(y 2 + z 2 ) 
cosh 2 (as) 



(4.1) 



where x is along the direction of charge flow. The param- 
eters A and B determine the strength of the potential 
in the longitudinal and transverse direction, respectively, 
and 1/a measures the extension of the potential along the 
x direction. Here, we are interested to study transport on 



4 




could only study longitudinal transmission near unity 
since its validity is limited for V(x, y, z) -C E F . 

We choose the incident wave-function k (£,?/, z) with 
three components of initial momentum k x ,k 3 y , k 3 z as a 
combination of plane-waves with proper boundary con- 
ditions, i.e., open boundary condition along the direc- 
tion of transport and <pk(x,y = 0,L±,z) — <j)\^{x,y,z — 
0. /. : 0. 



4>k(x,y,z) = 



ik x x 



cos{k J y) cos(k 3 z z) , (4.3) 



FIG. 2: A schematic of the potential (|4.1[) - for a constant 
value of z - used in this work to represent the scattering off a 
nanoj unction. 



a length scale such that 1/a <C l ee ~ y/Dr^, with D be- 
ing the diffusion constant. This simple integrable model 
allows us to calculate both the elastic momentum relax- 
ation time T k , the longitudinal transmission coefficient 
r^ x (defined as the transmission probability for processes 
which do not change the incident transverse momentum) , 
and the equilibrium density n(x). In the linear- response 
regime, i\ due to boundary scattering from the poten- 
tial (|4.ip can be estimated by employing the generalized 
Fermi golden rule, 



— — L X L 2 



dk! 

4^h 



<5(e k ,-e k ) |(</> k |y(a^,z)|Vv)|f4-2) 



with L x ~ l ee along x and L± being the length of the elec- 
trodes in the direction perpendicular to x. We use e k = 
h 2 k 2 /2m e with m e being the electron mass. The general- 
ized Fermi golden rule allows us to study the full range of 
longitudinal transmission for all values V{x,y,z), while 
in the Born approximation of the Fermi golden rule one 



We use the exact scattering state <fk' x (%) of the one- 
dimensional potential 1/ cosh 2 (ax) along the x— direction 
and plane waves with the same boundary conditions 
as the incident state in the transverse directions. For 
the transverse wave-functions, plane- waves are a bet- 
ter choice than the harmonic potential eigenstates in 
the asymptotic region, since in that region the curva- 
ture of the harmonic potential becomes infinitesimally 
small. Surely, the use of plane waves - instead of an ex- 
act scattering state |V> k ') of the full potential in Eq. (|4.1[) 
- would seem a limitation of our calculations. How- 
ever, we stress that the use of an asymptotic state of 
^k'(x,y,z) is also not the right choice in Eq. (14.21) . 
since the asymptotic state is a better approximation 
only far away from the action of the potential (where 
V(x,y,z)=0), while \(^ k \V(x,y, z)\^)\ of Eq. @3) is 
nonzero only for V(x,y, z) ^ 0. 

The scattered wave-function ^ k '(x, y, z) with momen- 



tum k' x ,k"P, k v z is thus 



tpk>(x,y,z) 



ML 



-ip K {x)cos(kPy)cos(k p z z) , (4.4) 



with 



J 



<PK (*) = (1 - er ik ' j2a FHK/u - s, -ik'Ja + 8 + 1, -ik'Ja + 1,(1- 0/2] , 



(4.5) 



r 



and £ = tanh(ax), where 
kP' = (2p' + l)7r/L ± with p,p' 



kP y = 



(2p + l)7r/Lj_ and 
0,±1,±2.., s = (-1 + 



\fl — 8m e A/a 2 h 2 )/2 and F(j3,j,S, z) is the hypergeo- 
metric function. The normalization factor M is fixed by 



fX%dx\ lpk ,(x)\ 2 =Af. 



-ix/2 



We separately evaluate r k , r+ the longitudinal and 
transverse relaxation times, respectively, with l/r k = 

1/ T k + ^/ T )t ■ We include in r k ' only those contributions 
of momentum relaxation which do not change the index 
of incident transverse channels, while all other processes 
are included in r k . Due to its symmetry in the trans- 
verse plane, the potential in Eq. (|4.1I) does not allow 



processes in which the y and z components of the in- 
cident momentum change simultaneously. For example, 
if we choose the incident state with j, j' = 0, the allowed 
contributions in t£ come from the channels with indices 
P = 1, 2, • • • ,p c , p' = and p = 0, p' = 1, 2, • • ■ ,p c , 
where p c is the maximum number of allowed transverse 
channels. 

On physical grounds of local thermodynamic equilib- 
rium in the electrodes, we can assume that the incident 
state is prepared with the lowest transverse momentum 
component, i.e., = 0. Therefore, for an incident 
state at the Fermi energy, k x — kp^l — 2(Tr/k F L±) 2 
and ky — k° z — n/L±. The longitudinal transmission co- 
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FIG. 3: Longitudinal momentum relaxation time r k ' as a func- 
tion of the transmission coefficient Tk m along the a:-direction 
for r s = 3. Inset shows the corresponding transverse momen- 
tum relaxation time T k . 



efficient T^ x along x for the longitudinal momentum k x of 
the incident state is calculated from iff. (x), and is given 
bylH 



sinh (jrk x /a) 



(4.6) 



sinh 2 (7rfc x /a) + cosh 2 (§ v / 8m e A/h 2 a 2 - 1) 
with 8m e A/a 2 h 2 > 1. We find after some algebra 

' "'■ ■l(-k^(A + BLi^)\ (4.7) 



h 3 \t , 

m e B 2 L 4 



' 6tt 2 
I (2p+l) 2 



^ 3 s=± iV=i,2.. |fc " |p4(p+1) 

where J(fcP) = (I/a/AT) /^f/ 2 dxeT ik * x ip k v (x)/ cosh 2 (ax). 

fcP [= fc Fv /l - ((2p + 1) 2 + l)7r 2 /(fc F L_L) 2 ] is the lon- 
gitudinal momentum of the scattered electron in the 
transverse channel of indices p and p' = 0. 

The transverse and longitudinal relaxation times are 
plotted in Fig. [3] for the electron density r s — 3, and 
typical dimensions of nanoscale junctions, namely L x = 
20 nm, L± = 5 nm and 1/a = 2 run. We always keep 
the transverse confinement potential strength (BL\ > 
55 eV) at far above the Fermi energy (e_p = 5.57 eV 
for r s = 3), and tune the height (A) of the longitudinal 
barrier across the Fermi energy to have different values 
of the transmission coefficient . 

Figure [3] confirms our initial hypothesis, namely that 
the elastic relaxation time is the dominant contribution 
to the total relaxation time for |k| = kp and for a wide 
range of temperatures: f/r = l/r ee + 1/Vk ~ We 
finally note that if we start with an incoming state with 
i> j' 7^ 0j then the incident longitudinal momentum k x at 
the Fermi energy is smaller than that of the incident state 
j, j' = we have considered so far. The corresponding 
longitudinal transmission Tfe x is also relatively smaller for 
the same longitudinal barrier height A. Thus, we find a 



FIG. 4: Zero- frequency xc viscosity r] xc in a nanojunction as 
a function of the transmission coefficient for r s = 3. 



faster longitudinal relaxation and a slower transverse re- 
laxation time at the same barrier height. However, if we 
plot the total relaxation time as a function of the longi- 
tudinal transmission coefficient by reducing the barrier 
height, we find that the total relaxation time is qualita- 
tively and quantitatively similar as for the — initial 
state. 

We are now ready to estimate the desired viscosity 
from Eq. (|3.2p . and from Eq. (|2 . 1[) the dynamical correc- 
tions to the single-particle resistance. For this we need 
only the xc component of the shear modulus, since the 
xc viscosity in the d.c. limit is r) xc (0) — (i xc ,<x>tf- This 
quantity is plotted in Fig. [4] with (i xc ,oo given by the 
second term on the right hand side of Eq. (|3.3|) and tf 
evaluated above for different values of Tk x . 

Finally, since in the present model the transverse vari- 
ation of the density is small, we evaluate only the lon- 
gitudinal contribution to Rdyn from Eq. (|2.I|) using the 
full wavefunction ^(x^y, z) and by defining the planar 
average 



(n(x)) = jy 



L ± /2 r L ± /2 



X J-L ± /2J-L ± /2 



dydz \Mx,y, z)\l4.9) 



k<k f 



with the proper normalization of (n(x)) corresponding 
to the bulk density of r s = 3 deep into the leads. This 
quantity is plotted in Fig. [5j Introduction of the poten- 
tial V{x,y,z) induces two effects: one is the reduction 
of the xc viscosity entering Eq . (|2.I[) due to fast momen- 
tum relaxation, and the other is the increase of the den- 
sity gradient. These two effects together determine the 
value of the dynamical resistance corrections in Fig. [5] 
In agreement with what was previously reported [7j , this 
dynamical resistance is relatively small at large trans- 
missions due to the small variation of the density across 
the junction. However, it increases substantially at very 
low transmissions with values that can greatly exceed, 
for this particular model, the dynamical resistances esti- 
mated previously [f~7T |. 
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FIG. 5: Dynamical viscous resistance as a function of trans- 
mission coefficient for the potential (|4.f [I , The inset shows 
the averaged density {n(x)) across the nano-junction for the 
longitudinal transmission coefficient 0.1 and 0.9. Other pa- 
rameters are in the text. 

V. CONCLUSION 

In conclusion, we have introduced a dispersion-relation 
approach to estimate the d.c. viscosity which provides 



very good agreement with the estimates obtained us- 
ing the standard non-perturbative calculation for the 
homogeneous electron liquid. This approach tremen- 
dously simplifies the calculations of the viscosity and 
allows to estimate this quantity in the presence of a 
nano-constriction where the momentum relaxation time 
is dominated by the elastic collisions at the junction. 
We have then computed the many-body contribution to 
the resistance of the junction for an integrable potential 
and found that while this resistance is relatively small 
for transparent barriers it is substantially higher for low 
transmission barriers, a fact which goes in the right di- 
rection in explaining the well-known (and yet unsolved) 
discrepancy between theory and experiments in molecu- 
lar junctions [r3 |. 
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